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EXPLOSION PHENOMENA IN STOCHASTIC 
COAGULATION FRAGMENTATION MODELS 

By Wolfgang Wagner 

Weierstrass Institute for Applied Analysis and Stochastics 

First we establish explosion criteria for jump processes with an 
arbitrary locally compact separable metric state space. Then these re- 
sults are applied to two stochastic coagulation-fragmentation models — 
the direct simulation model and the mass flow model. In the pure co- 
agulation case, there is almost sure explosion in the mass flow model 
for arbitrary homogeneous coagulation kernels with exponent bigger 
than 1. In the case of pure multiple fragmentation with a continu- 
ous size space, explosion occurs in both models provided the total 
fragmentation rate grows sufficiently fast at zero. However, an ex- 
ample shows that the explosion properties of both models are not 
equivalent. 



1. Introduction. Coagulation-fragmentation models are used in differ- 
ent application fields ranging from chemical engineering (reacting polymers, 
soot formation) or aerosol technology to astrophysics (formation of stars and 
planets). These models describe the behavior of a system of particles that 
are characterized by their sizes and move in a certain medium. The size of 
a particle changes either by coagulation (merging with another particle) or 
by fragmentation (breakage, splitting into pieces). We refer to the survey 
paper [1] for more details and references. 

Deterministic coagulation-fragmentation models are nonlinear evolution 
equations governing the macroscopic behavior of the particle system. In the 
case of a discrete size variable and spatial homogeneity, the coagulation 
equation takes the form 
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DC 

y=l 

where t > and x = 1, 2, Equation (1.1) goes back to Smoluchowski [29] 

and describes the time evolution of the particle number density (relative 
number of particles of a given size) . The coagulation kernel K is determined 
by the physics of the driving medium. It was observed that a certain kind 
of phase transition occurs when K grows sufficiently fast in its arguments. 
This phenomenon is called gelation and corresponds to a loss of mass in the 
macroscopic equation. The gelation point is defined as 

oo 

(1.2) tgd = inf{i > : mi(t) < mi(0)} where mi (t) = xcjt, x) . 

x=l 

Gelation is interpreted as the formation of infinitely large clusters in finite 
time. The presence of such clusters is not reflected in (1.1). 

Stochastic coagulation models go back to [14, 20, 22]. They are based on 
systems of particles 

(1.3) xi n \t), i = l,...,N^(t), N^(0)=n, 

which coagulate according to appropriate rates determined by the kernel K. 
These systems approximate the solution of the coagulation equation (1.1) 
in the sense 

oo 1 JVW(t) 

(1.4) 5>(x)c(t,s)= lim - E <p(*i n Ht)), t>0, 

* — ' n — >oo fi * — » 

x=l i=l 

for appropriate test functions (p (sequences in the discrete case). We re- 
fer to [8, 16, 26, 27] concerning rigorous results. An alternative stochastic 
coagulation model is related to the mass flow equation 

—c(t, x) = V — c(t, x - y)c(t, y) 

ot y 

y=l & 

(1-5) 

_^K( M ) 

where c(t,x) =xc(t,x). Equation (1.5) describes the time evolution of the 
particle mass density (relative mass of particles of a given size). A corre- 
sponding mass flow process is represented by a system of particles 

(1.6) x\ n \t), i = l,...,N<- n \t), N^(0)=n, 

with an appropriately modified evolution rule so that the solution to equa- 
tion (1.5) is approximated. Due to the equivalence of the equations, the 
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mass flow model provides an alternative approximation to the solution of 
the original coagulation equation (1.1), namely, 



(1.7) 



oo oo (p(x\ 

f(x)c(t, x) = J2 £ (i, x) 

2=1 X = l 



: lim - Y t)} , t>0. 

n^oo n f-f ~{n) 
i=l ■i'i 



We refer to [9] concerning a rigorous proof. A time-discrete version of the 
mass flow process was introduced and studied in [2]. 

In the direct simulation model (1.3), the mass of the system is conserved 
so that (1.4) implies 



oo , iV(")(t) n 

J2xc(t,x)~L £ *5»>(t) = ±5>W(0), t>0. 



n f— : n . 

X = l 1=1 1 = 1 



Thus, there is no way to approximate the gelation point following the defi- 
nition (1.2). Instead, the gelation effect is related to the formation of a very 
big cluster (order of the whole system) in finite time. In the mass flow model 
(1.6), the situation is different. According to (1.7), one obtains 

Vxc(i,x t>0, 

~ n 

x=l 

that is, the mass of the solution to equation (1.1) is approximated by the 
normalized number of particles in the mass flow system. Since the particles 
grow, there is a chance that their number drops due to some explosion 
phenomenon (infinitely many jumps in finite time). It was conjectured in [9] 
that, in case of gelling kernels, the mass flow process explodes and that the 
gelation time (1.2) is the limit (as n — > oo) of the (random) explosion times 
of the approximating finite particle systems. 

The paper contains three main results. First some rather general ex- 
plosion criteria for pure jump processes are obtained. Using ideas from 
[18, 19], previously known results for one-dimensional processes are gen- 
eralized to processes with a locally compact state space. This allows one to 
treat many stochastic coagulation-fragmentation models. Second we prove 
that the mass flow process explodes almost surely for a very wide class of 
gelling kernels in the case of pure coagulation. This result confirms the first 
part of the conjecture mentioned above. Its simple proof illustrates the use- 
fulness of the general explosion criteria. Third we study the case of pure 
multiple fragmentation for a continuous size space. Applying the general 
criteria, we prove explosion results for both the direct simulation model and 
the mass flow model, when the fragmentation rate at zero grows sufficiently 
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fast. These results correspond to another kind of phase transition — a loss of 
mass to zero (while gelation corresponds to a loss of mass to infinity). This 
phenomenon was studied a long time ago in [12] and named "transformation 
into dust." Later it was called "shattering" in the physical literature [23, 30]. 

The paper is organized as follows. In Section 2 we introduce the mini- 
mal jump process with an arbitrary locally compact separable metric state 
space and prove several sufficient criteria for explosion of this process. We 
give some examples and reproduce the previously known conditions for the 
one-dimensional case. In Section 3 we define the direct simulation model 
and the mass flow model. These models are extended to include both co- 
agulation and multiple fragmentation, as well as source and efflux terms, 
which are important in many applications. The relation of the stochastic 
models to the corresponding macroscopic equations is sketched heuristically, 
referring to convergence results available in the literature. In Section 4 we 
apply the explosion criteria from Section 2 to the stochastic coagulation- 
fragmentation models introduced in Section 3. First we consider the direct 
simulation model and prove a simple nonexplosion criterion in the general 
case, as well as an explosion result in the pure fragmentation case. Then we 
consider the mass flow model. In the pure coagulation case, we prove that 
there is almost sure explosion for arbitrary homogeneous coagulation ker- 
nels with exponent bigger than 1. Finally, we derive an explosion result in 
the pure fragmentation case and show that the explosion properties of both 
models are not equivalent. Section 5 contains some concluding remarks. 

2. Explosion criteria for jump processes. 

2.1. The minimal jump process. Let E and E' be separable metric spaces. 
The sets of measurable and continuous functions on E are denoted by 
M(E) and C(E), respectively. Furthermore, Mb(E) and V(E) are the sets of 
bounded Borel measures and of probability measures on the Borel-cr-algebra 
B{E). Finally, let 1b denote the indicator function of a set B, and 5^ be 
the Dirac measure on £ £ E. A kernel from E to E' (on E if E = E') is a 
function q:Ex B(E') — > [0,oo) such that 

q{-,B)€M(E) VBG B(E') and q{£, •) G M b (E') V£ € E. 

A kernel q is called compactly bounded if 

supg(£, E') < oo for any compact C C E. 

Sec 

Let q be a compactly bounded kernel on a locally compact separable 
metric space E. Let Coj Ci> • • • be a Markov chain in E with initial distribution 
uq £ V(E) and transition function p: E x B{E) — ► [0, 1] defined by 

(2-1) p(t,B) = l A(0 - A ^ j>U ' 

UM£):A(£) = 0, 
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where 

(2.2) A(£) = 

Let Tq,T\,... be independent and exponentially distributed random vari- 
ables with mean 1 that are also independent of (Cfc)> all defined on some 
probability space (Q,J-,F). Introduce the jump times 

11 T k 

t = 0, r « = E\77~T' 1 = 1,2,..., 
where Xfc/0 := oo, and the explosion time 
(2-3) ^ = Um 77 = ^— A_ . 

The minimal jump process, corresponding to the kernel q and the initial 
distribution vq, is defined as (cf. [11], page 263 and [25], page 69) 

(2.4) ^)-{Wtl™' 

where A ^ E determines the one-point compactification of E (cf. [3], page 205). 
The process is called regular, if 

(2.5) P(too = oo) = 1. 

Otherwise, the process is called explosive. Concerning the history of the 
subject we refer to [6], Prologue. 

Regularity and explosion of the minimal jump process can be studied 
using the properties (cf. [25], page 71) 

OO / OO rp \ 

(2.6) E-<°° => P(E— <°° =1 



fc=o° fc \k=o ak 



and 



oo 



(2.7) E^ = °° => P (E?=°° )=i. 



for any nonnegative sequence (afc). The independence of (7*.) and (Cfc) allows 
one to conclude from (2.6) and (2.7) that [cf. (2.3)] 



< OO 



°° I \ I I °° 1 \ 



which implies 



OO 



{tqo < oo} = < E 777 < oo > almost surely. 
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In particular, a necessary and sufficient condition for regularity (2.5) is (cf. 
[5], page 337) 

Correspondingly, the process is explosive if and only if 

(, 8 , <t^H>°- 

Note that boundedness of the waiting time parameter A implies regularity, 
and that 

OO -y 



(2.9) £ <00 \im( k = A. 

Thus, a necessary condition for explosion is 



(2.10) P lim Cfe = A >0. 

\k— >oo J 

2.2. Explosion criteria. First we formulate a general result concerning 
explosion of the minimal jump process defined in the previous section. 

Theorem 2.1. Let q be a compactly bounded kernel on a locally compact 
separable metric space E and (£&) be the corresponding Markov chain [cf. 
(2.1) and (2.2)]. Consider the sets 

(2.11) E £ ( V ) = ^eE: J E [v(£i)-Ti(0]q(S,d£x)>e} 
and 

(2.12) n £ (7]) = {ujen:Ck(u)£Es(v)i Vk>k(u), for some k{u)}, 
where r] is a bounded measurable function on E and e > 0. Then 

(2.13) 12^77^ < 00 a.s. onn £ (r]), Ve>0, 

that is, the minimal jump process, corresponding to the kernel q, explodes 
almost surely on the set of all trajectories staying in E e {rj) for sufficiently 
large k and some e > 0. 

Remark 2.2. According to Theorem 2.1, a sufficient condition for ex- 
plosion is the existence of some bounded measurable function r/ such that 

(2.14) F(Q S (7))) > for some e > 0. 
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It follows from (2.9) and (2.13) that 

lim = A a.s. on Q £ (rj), Ve > 0. 

k— >oo 

Thus, condition (2.14) puts some restriction on the possible structure of the 
set (2.11). In particular, this set can not be contained in any compact. 

Before proving Theorem 2.1, we derive several corollaries providing suffi- 
cient conditions for explosion. 

Corollary 2.3. Suppose that 

(2.15) /[r?(6)-r?(0M^^i)>£ V£G£+ 

JE 

for some bounded measurable function r/, some subset E + C E and some 
e>0. Then 

OO ^ 

V — — T < oo a.s. on {Ck G E + , Vfe}, 
k=o A ^ 

that is, the minimal jump process explodes almost surely on the set of all 
trajectories living in E + . 

Proof. Assumption (2.15) implies E + C E e (rj) and {( k G E + , V/c} C 
O e (r/) so that the assertion follows from (2.13). □ 

A particularly simple criterion for explosion is obtained in the case E+ = 

E. 

Corollary 2.4. Suppose there exists some bounded measurable func- 
tion T] such that 

(2.16) [[v{Zi)-v{Z)MS,d£i)>e V^eE, 

JE 

for some e > 0. Then 



that is, the minimal jump process explodes almost surely, for any initial 
distribution. 

Remark 2.5. Under the assumptions of Corollary 2.4, it follows from 
(2.9) that 

'lim Cfc = a) = 1- 

k — >oo / 
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Corollary 2.6. Suppose there exists a bounded measurable function rj 
such that 

(2-17) / [Ti(Si)-ri(t)W,dZi)>e V££C, 

J E 

for some compact C C E and some e > 0. Then 

00 1 f 1 

(2.18) X!t7FT <00 a.s. on <Mim Cfc = A L 

^ A (Cfc) U-* 00 J 

i/iat is, i/ie minimal jump process explodes almost surely on the set of all 
trajectories tending to the compactification point. 

Proof. Assumption (2.17) implies E\CdE E (rj) and {lim^oo = 
A} C £le(rj) so that the assertion follows from (2.13). □ 

Remark 2.7. Under the assumptions of Corollary 2.6, it follows from 
(2.9) and (2.18) that 

{£^ <oo rfe° a=A } a& 

so that the necessary condition (2.10) is also sufficient for explosion of the 
minimal jump process. Note that trajectories may stay in the compact set 
C forever. 

The proof of Theorem 2.1 is prepared by the following lemma. 
Lemma 2.8. Let n be a bounded measurable function on E. Then 

oo 

3 ^[E(?7(Cfc+i)|Cfc) -r)(Ck)] fi^ a.s. on O ("), 

k=0 

that is, the infinite sum is finite for almost all oj G ilo( r ?) [c/- (2.12)]. 
Proof. The sequence 

n-l 

(2.19) W n = J2mv(Ck+i)\tk)-v((k)]-ri((n), n>l,W = -r)(( ), 

fe=0 

is a martingale with respect to the filtration of (Cfc)- Representing it in the 
form 

W n = J2 mv(Sk+l)\(k)-ri((k)] 

k=0,...,n-l:C k eE (v) 

(2.20) 

+ [Hri((k+I)\tk)-V(tk)}-V(tn) 
k=0,...,n-l:Ck<£E (r)) 
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and introducing the sequence of stopping times 

a N = M{k > div-i : Cfe i E {v)}, N>2, 
CTl = inf{^>0:a^^o(r/)}, 
one concludes that 
W min(n , aN) > -sup|77(OI[2#{fc = 0, . . .,mm(n,a N ) - 1 : Q k i E Q (rj)} + 1] 

> — (2JV — 1) sup |»7(0I Vn,JV>l. 

Since the stopped process W min ^ n ^ N -j is a martingale and bounded from 
below, it has a.s. finite limits (cf., e.g., [28], Sections 11.49 and 57), that is, 

3 lim W min(n (JJV ) finite ) =1 ViV. 

Consequently, one obtains 

(2.21) 3 lim W n finite a.s. on \a N = oo} ViV. 

n— >oo 

Note that the first sum at the right-hand side of (2.20) increases, while the 
second sum has at most N — 1 elements (if cr/v = oo) and the last term is 
bounded. Thus, (2.21) implies 

(2.22) 3 lim rj(( n ) finite a.s. on {aN = 00} ViV. 

n — >oo 

Finally, one concludes from (2.19), (2.21) and (2.22) that 

00 

3 ^[E(r/(Cfe + i)|Cfe) ~v(Ck)] finite a.s. on {a N = 00} ViV, 

k=0 

and the assertion follows from the fact that Flo(r]) = {JnWn = 00}. □ 
Proof of Theorem 2.1. Note that 

[V(ti) ~ T)(0]q(Z,dti) = A(OpEMCfc +1 )|a = " 



and 



A(Cfc) " e 

Since 

fe-l 



< ipE(7/(C fc +i)|Cfc) - v(Ck)} if C* e EM- 



LXJ -1 K — 1-1 -1 LXJ 

Et77T^Et77T + -E PE(r?(Cfc + i)|G) - 17(a)] on n^fa), 
the assertion follows from Lemma 2.8. □ 
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2.3. Examples. 

Example 2.9. Consider the case 

£ = {1,2,...}, g(£,dei) = A(0%u(#i). Co = l- 

The trajectory of the Markov chain is deterministic, ^ = k + 1, k > 0, and 
the necessary and sufficient condition for explosion (2.8) takes the form 

(2.23) g^<oo. 

First we show that the sufficient condition (2.16) of Corollary 2.4 can always 
be satisfied. Indeed, choosing the function 

C-1 1 

^) = E;^)' ^ 2 < ^ 1 ) = ' 

which is bounded if (2.23) holds, one obtains 

fo(£ + l)-»7(0]A(0 = l V ^ Gi? - 
Furthermore, for A(£) = £ and any strictly increasing bounded function rj, 
this example illustrates that the condition 

/ [v(£i)-v{ZMS,dZi)>o vee^, 

JE 

instead of (2.16), would not be sufficient. Finally, for ??(£) =£ and constant 
A(£), this example illustrates that condition (2.16) for unbounded rj would 
not be sufficient. 

Example 2.10. Consider the one-dimensional case E= [1, oo) and the 
bounded measurable function 

v(0 = -£,~ a i £ G -E 1 for some a > 0. 
Condition (2.17) takes the form 

(2.24) A(f)f ~ q e(i - (^j | Co = t) > e for sufficiently large f . 
With the notation 

n„ Q (£)=E(\-(£) a |Co = £) and /(O^f* 1 , 

condition (2.24) transforms into 

A(£)£ra_a,(£) > /(£) for sufficiently large £, 
so that the result of [18], Theorem 3 is basically reproduced by Corollary 2.6. 
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3. Deterministic and stochastic coagulation fragmentation models. In 

this section we introduce two stochastic coagulation-fragmentation models 
that are covered by the framework of Section 2. Explosion phenomena in 
these models will be studied in Section 4. 

The two models are related to each other in the sense of (1.4) and (1.7). In 
order to demonstrate this, we let them depend on a parameter n = 1, 2, . . . 
and illustrate heuristically their relation (when n — > oo) to corresponding 
macroscopic equations of the type (1.1) and (1.5). Rigorous convergence 
results and more references can be found, for example, in [10]. The treatment 
of the macroscopic equations in the general setting (including coagulation, 
multiple fragmentation, source and efflux terms) is of independent interest, 
since many equations known from the literature (weak and strong, discrete 
and continuous) are covered in a unified way. 

We consider stochastic processes with the state space 

(3.1) EW = f±-Y;6 Xi :N>0, %i & X, i = 1,...,n\, 
where X is a locally compact separable metric space. A state £ € 

E {n) 

IS 

interpreted as a system of particles with types from X and weights ^ . Stan- 
dard examples are 

(3.2) X = {1,2,...} or X = (0,oo), 

corresponding to discrete and continuous particles sizes, respectively. How- 
ever, many of the constructions of this section are valid in the general setup, 
which allows one to cover particles consisting of several species, and other 
cases of practical interest. 

Remark 3.1. The space E^ n \ endowed with an appropriate metric of 
weak convergence, is locally compact and separable (cf. Remark 2.10 and 
Lemma 5.1 in [10]). Note that — > A in E^ if either — > oo or Xk,i — ► Ax 
for some i, where Ax denotes the compactification point of X. 

3.1. The direct simulation model. Consider the kernel 
qW{Z,B) 

= n f l B (J s (£,x))S(dx) 
Jx 

(3.3) N N 

+ £l B (J e (C,i))e(xi) + Y, / lB(MZ,i,z))F(xi,dz) 
i=i i=i Jz 

l<i^j<N 
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with the jump transformations 



(3.4) 



MM 


=£ + 


1, 

— d x , 

n 




= e- 


1, 

n 


J F (£,i,z) 


=£- 


n n k 






1 1 

-fc + 5 X ] + -4<+xj- 

n n 



Here 5 G .M&(X) is some source measure. The nonnegative function e G 
C{X) denotes the efflux intensity. The compactly bounded fragmentation 
kernel F from X to 



(3.5) 



Z= \Jx k 



k=2 



is assumed to satisfy the mass conservation property 

(3.6) F(x,Z\Z(x)) = Vi£l, 

where Z{x) = {z G Z : z\ + • • • + z\. = x}. The nonnegative function K G 
C{X x X) denotes the coagulation kernel. Under the above assumptions, the 
kernel (3.3) is compactly bounded (cf. [10], Lemma 5.1). The corresponding 
minimal jump process is called direct simulation process. 
Note that 

(3.7) =n S(X)+ [ e(x)£(dx) 

Jx 

+ [ F(x,Z)£(dx) + -^ E K ( x U x j) 



and 



E( n ) 



x 



(3.8) 



A 



+ 



l<ij*=j<N 

V (x)^(dx)- I ^(x)e(dx)lg( n )(e,dei) 
Jx 

ip(x)S(dx) — / ip(x)e(x)^(dx) 
Jx 

[tp(zi) H h p(z k ) - tp(x)]F(x, dz)i(dx) 



x Jz 
1 



+ 2^2 ZZ [f{Xi + Xj)-ip{Xi)-lf{Xj)}K{x i ,X j ), 



l<iy£j<N 
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for any £ G E 1 ^) and appropriate test functions <*p. Dynkin's formula and 

(3.8) suggest that the direct simulation process corresponds (for n — ► oo) to 
the macroscopic coagulation-fragmentation equation 

^| p(x)/i(i,da;) 

(3.9) +/ [ [<p(z 1 ) + --- + ip(z k )-<p(x)]F(x,dz)n{t,dx) 

Jx Jz 



+ o / / [y(a? + y) - <p(z) - <p(y)} 
* Jx Jx 

x K(x,y)n(t,dx)/j,(t,dy), t>0, 



with the initial condition fi(0,dx) = fio(dx), where //o is the limit of the 
initial state of the process. 

Let denote the restrictions of the fragmentation kernel F to the 

sets X k [cf. (3.5)]. Define the symmetrized fragmentation kernel F sym by its 
corresponding restrictions 

^sym(^' dz\ , . . . , dZfc) 

= jj J2 F{k) ( x ' dz Ai) ' • • • ' dz Hk) )' A; = 2, 3, ... , 

7T 

where the sum is taken over all permutations of {1, ...,k}. Introduce the 
1-marginals of the kernels (3.10) 

(3.11) fff ( x , dy) = Fgl(x, dy,X,..., X) 
and the kernel 

oo 

(3.12) F£l(x,dy) = Y,kF^(x,dy). 

k=2 

Then one obtains 

[tp(zi) H \-<p(z k )]F(x,dz) 



oo „ 

(3.13) = £ / Mz l ) + --- + V (z k )}F^l(x,dz l ,...,dz k ) 

k=2 JX 

oo 

= 5>/ <p{y)F^{x,dy)= <p(y)F^(x,dy). 



k=2 



X JX 
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Using the symmetrized coagulation kernel 

(3.14) K sym (x,y) = \[K(x,y) + K(y,x)}, 
one obtains 

\\ \ [^p(x + y)-^(x)-(p(y)]K(x,y)n(t,dx)fj,(t,dy) 
Jx Jx 

(3.15) = \ J J [y{x + y) ~ ~ ( p(y))K sym {x,y)fj,(t,dx)fi(t,dy) 

/ [\vi x + y)~ (p(x)]K sym (x, y)n(t, dx)n(t, dy). 
xJx 

Remark 3.2. According to (3.13) and (3.15), the macroscopic equa- 
tion (3.9) does not change if the fragmentation and coagulation kernels of 
the direct simulation process are replaced by their symmetrizations. 

Preparing the transition from (3.9) to an equation for the densities 

n(t, dx) = c(t, x) dx, 
we assume in the continuous case X = (0, oo) that 

(3.16) S(dx)=s(x)dx 
and [cf. (3.12)] 

(3.17) Fgl(x, dy) = /g,(i, y) dy. 
Note the identity 



(3.18) / / i[;(x,y)dydx= / / ip(x - y,y) dy dx, 

Jo Jo Jo Jo 

where ip is an appropriate test function. Using (3.13), (3.15) and (3.18), one 

obtains the continuous coagulation-fragmentation equation with source and 

efflux terms 

d 

-c(t, x) = s(x) — e(x)c(t, x) 



dt 





+ / fsl ) m ( x + y, x ) c ( t , x + y)dy -c(t,x)F(x,z) 



(3.1 9 ) 

+ 77 / K sym (x - y, y)c{t, x - y)c(t, y) dy 
z> JO 

•oo 

K syra [x, y)c(t, x)c(t, y) dy. 

In the discrete case X = {1,2,...}, an equation analogous to (3.19) is ob- 
tained, with integrals replaced by sums. 
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Remark 3.3. According to (3.13), the kernel (3.12) can be expressed 
in the form 

(3.20) F^l(x, dy) = f [5 Z1 (dy) + ■ ■ ■ + S Zk (dy))F(x, dz). 

Thus, the quantity 



^-^Fgl(x, (0, y)) = E x £ 1(0,,) W 



K x J2 1 

i-.Zi<y 

represents the average number of fragments with size less than y, resulting 
from the fragmentation of a particle of size x. In particular, the average 
number of fragments is expressed as 



(3.21) 



F(x,Z) F ^ m( - X,X) - 



Note that F^}a(x,[e,x)) < oo,Ve G (0,x), due to mass conservation (3.6). 
Thus, for test functions tp with compact support, the corresponding term in 
(3.9) is finite even if the average number of fragments is infinite. Finally, we 
mention that mass conservation implies 

F ( - k \x,X k )= [ (z 1 + --- + z k )F^(x,dz 1 ,...,dz k ) = k ( yF^\x,dy) 



XI 

and 
(3.22) 



F(x,Z) = - [ yFgl(x,dy). 
x j x 



Example 3.4. Consider the binary fragmentation case 

(3.23) F(x,dz) = F ( - 2) (x,dz 1 ,dz 2 ) = \F {1) (x,dz 1 )5 x - zi (dz 2 ), 
where F^ is a kernel on X such that 

F {1) (x,X\(0,x)) = VxGX 
One obtains [cf. (3.10)-(3.12)] 

Fgl(x,d Zl ,dz 2 ) = l[F (1 \x,dz 1 )5 x - Zl (dz 2 ) + F ( ' 1) (x,dz2)5 x - Z2 (dz 1 )] 
so that 

Fg^(x,dy) = l\F^(x,dy)+ f F« (x,dx 1 )S x ^ xl (dy) 

and 

(3.24) F^l(x,dy) = \\f^{x , dy) + [f^ (x,d Xl )5 x . xl (dy) 



x 
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Assuming 

F^(x,dy)=l^ x) (y)f^(x,y)dy, 
one obtains from (3.17) and (3.24) that 

f$Ux,y) = ±[fU(x,y) + fU(x,x-y)]. 

Since 



F{x,Z) = \ / fW(x,y)dy 



fsym(x,y)dy, 



the fragmentation term in (3.19) takes the usual form. Note that (3.22) 
implies 



1 -F^(x,X) = - I ' yFgl (x,dy) 



x Jx 



xeX. 



3.2. The mass flow model. The solution fi(t, dx) of equation (3.9) repre- 
sents the flow of concentration in the size space X. Since the total mass of 
the system is determined as f x x/j,(t,dx), we call the function 

jl(t,dx) = x/i(t,dx), t > 0, 

the mass flow. Considering test functions of the form (p(x) = xip(x), one 
obtains [cf. (3.9) and (3.13)] 

yP(y) F syb( x i d y)i J '( t ' dx )= / i>(y) F (x,dy)ii(t,dx), 



x Jx 
where 

(3.25) 

and [cf. (3.15)] 



x Jx 



F(x,dy) = ^l(x,dy) 



x Jx 



— {x + y)ip(x + y) — xtp(x) 



Ksym(x, y)n(t, dx)/j,(t, dy) 



x Jx 



[xtp(x + y)- x*li(x)]K sym (x, y)n(t, dx)n(t, dy) 



K sym (x,y) _ 



tp{x + y) -ip(x)]'- 
lx Jx' y 

Since, according to (3.22), 

(3.26) F(x,X) = F(x,Z), 



fi(t,dx)fi(t,dy). 
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equation (3.9) takes the form of the mass flow equation 



(3.27) 



x 



tp(x)xS(dx) 



i/j(x)e(x)fl(t,dx) +/ / [ip(y) - ip(x)]F(x,dy)fl(t,dx) 
x Jx Jx 



+ [ f Mx + y)- Ksym{X ' V) flit, dyMt, dx) , t>0, 
JxJx y 

for appropriate test functions ip and some initial condition. 

In the continuous case X = (0, oo), we assume (3.16) and (3.17). It follows 
from (3.25) that F(x,dy) = f(x,y)dy, where 

f(x,y) = ~,fSL{x,y). 

Using (3.18), one obtains from (3.27) an equation for the densities 

jl(t, dx) = c(t, x) dx, 
namely, the continuous mass flow equation with source and efflux terms 
d 

-c{t, x) = xs(x) — e(x)c(t, x) 



dt 



K sym (x V>v) 5{t )5{t )d 

y 



K sym (x,y) 



c(t,x)c(t,y)dy 



}{x + y,x)c(t,x + y)dy - / f{x,y)c{t,x)dy 



This equation is equivalent to (3.19). 
We assume 

(3.28) / xS(dx) < oo 

Jx 

and introduce the modified kernel [cf. (3.3)] 

B)=n f l B (J s (ti,x))xS(dx) 
Jx 

N N . 

(3.29) +J2^B{Je(Z,i))e{xi) + Y, / lB(MS,i>v))F(?i,<ht) 

i=l i=l Jx 

+ - f; t B (JK(Z,i,j)) Ksym{xi,Xj \ BeB( E (% 

i,j=l 
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with the modified jump transformations [cf. (3.4)] 
(3.30) 



Jf(Z, hV)=£- Sxi + -S y , 
n n 



n n J 

The minimal jump process, corresponding to the kernel (3.29), is called mass 
flow process. 

Following (3.20), the mass flow fragmentation kernel (3.25) can be repre- 
sented in the form 



(3.31) 



F(x,dy) 



x x 



F{x,dz). 



Since the total fragmentation rate does not change [cf. (3.26)], the "frag- 
ment" (next state) in the mass flow process is chosen as follows: first frag- 
ments z\ , . . . , Zk are generated according to the normalized direct simulation 
fragmentation kernel F(x, dz); then one of them is chosen with probabilities 
proportional to their individual masses. 

Note that [cf. (3.7), (3.8) and (3.28)-(3.30)] 



(3.32) 



and 



n 



x 



xS(dx) + / e(x)£ t (dx 
Jx 



x 



F{x,X)i{dx) + 



K sym (x,y) 



x Jx 



ip(x)^i(dx) — / ip(x)^(dx) 



x 



(3.33) 



x 



^(x)xS(dx) — / ip(x)e(x)^(dx) 



x 



+ 
+ 



X JX 



X JX 



[ifj(y)-ij(x)]F(x,dy)Z(dx) 

{x + y)-^x)] Ka ^ y \ {dx)i{dy), 



for any £ £ and appropriate test functions t/j. The corresponding macro- 
scopic equation (obtained for oo) is (3.27). 



3.3. Comments. The study of the coagulation-fragmentation equation 
(3.19) (with s = e = and binary fragmentation) goes back to [24]. 
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In the pure fragmentation case S = 0, e = 0, K = 0, the mass flow equa- 
tion (3.27) with test functions ip(y) = t^,x)(y) implies 



J^ft(0,s)) 



[l(o,z)(yi) - l-(o,x)(y)]F(y,dyi)(j,(t,dy) 



(3.34) 



x Jx 

/ F(y,dyi)fl(t,dy) 

x J0 

oo 

F(y,(0,x))d y /i(t,(0,y)). 



According to (3.31), one obtains 
F(y,(0,s)) = 



|^(0,*)) + --- + |^(M) 



F(y,Z)E, 



1 



E 



y ■ 

" 2 : Zi<x 



so that equation (3.34) is identical with equation (2) in [12]. Note that 
fl(t, (0,x)) is the average sum of masses of particles of mass smaller than x 
at time t. Since [cf. (3.26)] 



(3.35) 



E, 



7. E Z i 



y 



i : Zi<x 



F(y,(0,x)) 
F(y,X) ' 



the one-dimensional Markov process introduced in [12], page 279, is just the 
mass flow process. 

4. Explosion in stochastic coagulation fragmentation models. Here we 
apply the criteria from Section 2 to the models introduced in Section 3. The 
trajectories of the underlying Markov chain take values in the space (3.1). 
The state space of single particles is (3.2). The jump kernels and waiting time 
parameters are determined by (3.3), (3.7) and (3.29), (3.32), respectively. We 
fix the parameter n and skip the corresponding superscripts. 



Remark 4.1. Condition 

/ 00 

(4.1) P £ 



1 



implies 



' HCk) 



< 00 



< 00 



Co = e =1 V^eEcE, 



E 



< 00 



Co=eU(d£) = i, 



that is, explosion with probability one [cf. (2.8)], for any initial distribution 
uq on E. 
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It is often easier to prove (4.1), since trajectories starting at a given £ 
remain in a certain part of the state space, where some of the sufficient 
conditions can be checked. Note that, according to (2.10) and Remark 3.1, 
explosion implies that (with positive probability) either the number of par- 
ticles in the system reaches oo or some of the particles reach or oo. 

4.1. Explosion of the direct simulation process. We consider the process 
with the kernel (3.3). In the pure coagulation case the waiting time param- 
eter (3.7) takes the form 

Particle sizes increase so that they cannot reach zero. Both the growth of 
individual particles and the number of particles in the system are bounded 
due to mass conservation. So there is no explosion provided the coagulation 
kernel is bounded on compacts. Analogous arguments apply to the discrete 
coagulation-fragmentation case (without source term). Thus, all interest- 
ing cases in the sense of explosion should include continuous fragmentation. 
First, we study the situation, when the total fragmentation rate is bounded 
at zero, and provide some sufficient conditions for regularity. Finally, The- 
orem 4.3 gives a sufficient condition for explosion in the case of pure frag- 
mentation. 

Theorem 4.2. Consider the direct simulation kernel (3.3), where 

(4.2) F ik \x,X k ) =0 Vx£X,k>k F , 

(4.3) F(x,Z) <C F {l + x) VxGX, 

(4.4) K(x,y) <C K {x + y + xy) Vi,i/Gl, 

(4.5) e{x) <C e (l + x) Via, 
and 

(4.6) 5(Xn(C 5 ,oo)) = 0, 

for some constants k F , C F , Ck, C e and Cs- If there is no source term (Cs = 
0) or no coagulation term (Ck = 0), then the direct simulation process is 
regular. 

Proof. According to (4.3)-(4.6), the waiting time parameter (3.7) sat- 
isfies 



A(£) < nl S(X n (0, C s ]) + (C e + C F )[M (0 + Mi(£)] 

(4.7) 

■pMoCOM^O + MiCO 



2 
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where the notations: 

M ($)=Z(X), Mi(0= f xf(dx), ee£, 

are used. Consider trajectories (Co: Ci > • • • )• According to assumption (4.2), 
the number of particles in the system grows at most linearly, that is, 

(4.8) P(M (Cfc) < M (Co) + kk F , Vfc > 0) = 1. 

According to assumption (4.6), the mass of the system grows at most lin- 
early, that is, 

(4.9) P(Mi(Cfc) < Mi (Co) + kC s , Vfc> 0) = 1. 

If there is no coagulation term, then (4.7)-(4.9) imply (a.s.) 

, x A(C*) < n{S(X) + (C e + C F )[M (( ) + kk F + Mi (Co) + feC s ]} 
(4.10) 

Vfc>0. 

If there is no source term, then the mass of the system does not grow, that 
is, 

(4.11) P(Mi(Cfc) < Mi(Co), VA: > 0) = 1, 

and (4.7), (4.8), (4.11) imply (a.s.) 

A(Cfc) < n{{C e + G»[M (Co) + kk F + Mi (Co)] 



(4.12) 



+ ^[2[M (Co) + kk^M^Co) + Mx(Co) 2 ]} V/c > 0. 



In both cases (4.10) and (4.12), one obtains (a.s.) 



oo oo ^ 

so that regularity follows. □ 

Theorem 4.3. Consider the direct simulation kernel (3.3) in the pure 
fragmentation case (S = 0, e = 0, K = 0). If 

(4.13) F(x,Z)>^ \fx£X, for some C F > 0,a > 0, 

then the direct simulation process explodes almost surely, for any initial dis- 
tribution on the set E = {£ G E : N > 1}. 
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Proof. In order to make use of Corollary 2.3, we introduce the sets 

(4.14) E + (C) = UeE:-Yxi<C,N>l\, C > 0, 

I n T^i J 

and the function ry(£) = G E 1 , where 

f ( x ) = T~l — fl' x - °' < ^ - min (" 5 !)> 
1 + x p 

is an increasing bounded function. The derivative 
is decreasing so that 

(4.16) g(y) - g(x) > g \y){y - x) Vy>x>0. 

Furthermore, we note that Jensen's inequality implies (the convex function 
is x~ a , and the random variable is uniformly distributed over N different 
points) 

1 N ( 1 N \ ~ a 

(4.17) ^E^ r °> hvE^ V^>0,z = l,...,7V,iV = l,2,.... 

i=l \ i=l I 

According to (3.3), (3.4), (4.13) and (4.15)^(4.17), one obtains [cf. (4.14)] 

>(6)-??(0M£,^i) 

N 



(4.18) 



= E / Wf{^z)) -rj{0]F{^dz) 
i=\ Jz 

N 

>E / [9{{N+l)/n)-g{N/n)\F{x i ,dz) 

r N 
>-*V((JV+l)/n)$> 



i=l 

The order of N satisfies (3 — 1 — 2(3 + \ + a = a — {3>0, which makes the 
right-hand side of (4.18) bounded from below by some e > 0. Thus, condition 
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(2.15) is fulfilled and Corollary 2.3 implies almost sure explosion on the set 
of all trajectories (Co? Ci ; • • •) living in E + (C). Note that any £ G E satisfies 
£ G E + (C) for sufficiently large C, and 

Co£E + (C) =► Ck€#+(C) Vfca.s., 

since mass is conserved and the number of particles increases. Thus, the 
assertion follows from Remark 4.1. □ 



4.2. Explosion of the mass flow process. We consider the process with 
the kernel (3.29). In this model the growth of the size of individual particles 
is not bounded so that explosion is possible in the pure coagulation case. 
Theorem 4.4 gives a surprisingly easy and rather complete solution to this 
problem. As to pure fragmentation, again only the continuous case is of 
interest in the sense of explosion. Fragmentation in the mass flow process 
does not lead to a blow-up of the number of particles, but only to a decrease 
of their sizes. Theorem 4.7 provides sufficient conditions for explosion in the 
pure fragmentation case. 

Theorem 4.4. Consider the mass flow kernel (3.29) in the pure coagu- 
lation case (S = 0, e = 0, F = 0). Assume 

(4.19) K(x,y)>K(x,y), 
where K is homogeneous with exponent a > 1, that is, 

(4.20) K(cx,cy)=c a K(x,y) VoO.ijEl, 

and such that K(l, 1) > 0. Then the mass flow process explodes almost surely, 
for any initial distribution on the set E = {£ G E : N > 1}. 

Proof. In order to make use of Corollary 2.3, we introduce the sets 

E+(C, L) = jf G E : .jnm^i >C,N = l\, C > 0, L = 1, 2, . . . , 

and the bounded measurable function 

(4.21) V(0 = {L H{X) ^ if ^ E+ ^ L ^ UE, 

\ 0, otherwise, 

where 

H(x) = —x~P for some < /3 < a — 1. 
Consider trajectories (£o ; Ci> • • • )• Note that 

(4.22) ( eE + {C,L) ( k eE + (C,L) VA; a.s., 
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since the sizes of individual particles increase and the number of particles is 
conserved. According to (3.33), (3.14), (4.19) and (4.20), one obtains 

JE 

1 v— v r —8 I | \—0-i-Ksym{ x ii x j) 

= iLK x '+ x i) — 

n 2 .4-'/ sc,- 

1=1 

1 9-/3 L 

=^J-*(M)E*r 1HS 

i=l 

1 - 2~< 3 - 

> ^ i^(l,l)LC ,a ~ 1 ~ /3 >0 \/£EE + {C,L). 

Thus, condition (2.15) is fulfilled and Corollary 2.3 implies almost sure ex- 
plosion on the set of all trajectories living in E + (C,L). Note that any £ G 
satisfies £ G E + (C,L) for sufficiently small C and some L. Taking into ac- 
count (4.22), the assertion follows from Remark 4.1. □ 



Example 4.5. Consider the mass flow coagulation process MF(1) start- 
ing with one particle of size one, that is, C,q = 5\,n = 1. In this case the 
Markov chain is deterministic, 

Ck = $x k i x k = 2 , h = 0, 1, ... , 

and the sequence of waiting time parameters takes the form [cf. (3.32)] 

K(2 k ,2 k ) 



A(Cfc) 



2 k 



Assuming K = K and using (4.20), one obtains A(Cfc) = 2 k ^ a - 1 \ k = 0, 1, . . . , 
so that 

^- < co <^=^ a > 1. 

Example 4.6. Consider the mass flow coagulation process MF(2) start- 
ing with two particles of size one, that is, £o = 6i,n = 2. In this case the 
Markov chain takes the form 

C,k = \[5 Xk +S yh ], /c = 0,1,..., x = 2/0 = 1, 
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and the sequence of waiting time parameters is [cf. (3.32)] 



K(x k ,x k ) K(y k ,y k ) K(x k ,y k ) K(x k ,y k ) 



(4.23) A(C fc ) = ^ 

2 L x k y k y k x k 

If 

(4.24) K(x,y) > {xyf \/x,y£X, for some /?, 
then one obtains 

(4.25) A(G0 > + V r i + 4y^ +4~V k l 

We consider several special trajectories, for which the explosion property 
can be checked explicitly. 

First there is a trajectory of "fastest growth," 

(4.26) x k = 2 k , y k = l, k = 0,l,..., 

when each jump consists in doubling the first particle. If the coagulation 
kernel satisfies (4.24), then one obtains from (4.25) 

(4.27) A(C fc ) > \4 
and 

CO oo 2 

£ T77T ^ £ < °° if /3>0- 
fc =o A (CfcJ fc=0 l ' 

Thus, there is explosion on the trajectory (4.26) even for many nongelling 
kernels. 

Next we consider a trajectory of "second fastest growth," 

(4.28) (1, 1) - (1, 2) - (3, 2) - (3, 5) - (8, 5) -> (8, 13) - ■ • ■ , 

when self-interaction is avoided and alternatingly either the first particle is 
added to the second or vice versa. This trajectory is related to the Fibonacci 
numbers 



2 k y/5 



Integer 



V5\ 2 ) 



= 1, 2 



where Integer [a] denotes the nearest integer of the number a. Here both 
particles x k ,y k grow as C k for some C > 1. If the coagulation kernel satis- 
fies (4.24), then (4.25) implies \(Ck) ^ C^ 2 ^~^ k so that there is explosion on 
the trajectory (4.28) in the gelling case (3 > 1/2. 

Finally, we consider a trajectory of "slowest growth," 

(4.29) x k = l + k, y k = l, A; = 0,1,..., 
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when each jump consists in adding the second particle to the first one. If 
the coagulation kernel satisfies (4.24), then the estimate [cf. (4.27)] 

oo oo 2 

implies explosion in the case > 1. If, however, K{x,y) < xy, V x,y G X, 
then one obtains [cf. (4.23)] \(Ck) < %k + Vk = & + 2 and 

Thus, there is no explosion on the trajectory (4.29) even for many gelling 
kernels. 

Theorem 4.7. Consider the mass flow kernel (3.29) in the pure frag- 
mentation case (S = 0,e = 0,K = 0). Assume 

(4.30) F(x, X) > ^ V x £ X, for some C F > 0, a > 0, 
and 



(4.31) TT^R F(x,dy)< 7 <l WxeX. 
F(x,X) Jx \xj 

Then the mass flow process explodes almost surely, for any initial distribution 
on the set E = £ E : N > 1}. 

Proof. In order to make use of Corollary 2.3, we consider the sets 

E + (C,L) = UeE: max Xi<C,N = L\, C > 0, L = 1, 2, . . . , 
i=l,..., N J 

and the bounded measurable function (4.21) with H{x) = —x a . Consider 
trajectories (Co, Ci> • • • )■ Note that 

(4.32) Co eE + (C,L) => ( k £E + {C,L) Vfca.s., 

since the sizes of individual particles decrease and the number of particles 
is conserved. According to (3.33), (4.30) and (4.31), one obtains 

= -£ / K-y a ]F{x t ,dy) 



n 
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Thus, condition (2.15) is fulfilled and Corollary 2.3 implies almost sure ex- 
plosion on the set of all trajectories living in E + (C,L). Note that any £ G E 
satisfies £ G E + (C,L) for sufficiently large C and some L. Taking into ac- 
count (4.32), the assertion follows from Remark 4.1. □ 

In terms of the direct simulation fragmentation kernel F [cf. (3.25) and 
(3.26)], assumption (4.30) takes the form (4.13), while assumption (4.31) 
takes the form 

where the first factor represents the average number of fragments [cf. (3.21)]. 
The following examples illustrate assumption (4.33). 

Example 4.8. Consider the case of uniform binary fragmentation [cf. 
(3.23)] 

F^(x,dy) = F s « (x, dy) = F(x)U (0tX) (y) dy, x G X = (0, oo), 

where the function F determines the waiting time parameter. Assumption 
(4.33) is fulfilled, since 

2 r x „j.i . 2 



[ X y a + 1 dy = — -<1 Va>0. 
Jo a + 2 



Example 4.9. Consider the case of deterministic binary fragmentation 
[cf. (3.23)] 

(4.34) FW(x,dy) = F(x)5 K{x) (dy), 

where 

k(x)e(0,x) ViGI = (0,oo), 

and the function F determines the waiting time parameter. Note that [cf. 
(3.24)] 

FgUx,dy) = ?P- 

(4.35) 

F(x) 

= —j— [ s n(x)(dy) + $x-K(x)(dy)]- 

Thus, assumption (4.33) takes the form 

(4 .3 6 ) (^)" +1 + (^) <>+1 < 7<1 V * 6 X, 



8 K (x){dy)+ / $K(x)(dxi)5 x - Xl (dy) 



x 
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and is fulfilled, if 

71 x < k(x) < 72X V x G X, for some < 71 < 72 < 1, 

since 

K[x)_\ a+1 + / x-k(x) \ q+1 < «(s) a | z-/c(a;) 
x/ \ x / — x x 

<max(7 2 Q ,(l- 7l r)<l. 

In the last two examples we consider the case of deterministic binary 
fragmentation (4.34) and the mass flow fragmentation process with n = 1 
and one initial particle of size xq. Note that [cf. (3.25) and (4.35)] 

F(x) 

F(x,dy) = -^—[yS K ( x )(dy) + y5 x - K ( x ){dy)] 
F(x) 

= -2^-[K(x)S K (z)(dy) + [x- K(x)]5 x _ K{x) (dy)} 

and 

k(x)' 



-F(x,dy) = ^S K{x) (dy) + 



1 



8 X -K(x)(dy)- 



F(x,X) x 
Thus, the Markov chain takes the form ^ = 5 Xk , k = 0, 1, . . . , where 

K(x k ) 



(4.37) x k+l 



K(xfc), with probability 

Xk 

Xk — n(xk), with probability 1 



Xk 

The sequence of waiting time parameters is [cf. (3.32)] 

(4.38) \(( k ) = F(x k ,X) = lF(x k ). 

Example 4.10. Consider the case (4.34) with the function 

( x 1 ., 1 

- + -, if x> ± 

(4.39) k{x) = 1 2 4' 2 ' 

y — , otherwise. 

Note that assumption (4.33) [cf. (4.36)] is violated, since lim e ^o ^(5 + e ) = |- 
If xo > A , then the sequence 

2xq + 2 fc - 1 
% = 2^+1 ' fc = ° 5 1 '---. 



satisfies 



. 2x + 2 fc - 1 , 1 
K Wk) = ^+2 + 4 = 
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and corresponds to the trajectory of "slowest decrease" [cf. (4.37)]. There is 
no explosion on this trajectory, since it does not reach zero. However, it has 
a nonzero probability, namely, 



This example illustrates that some additional restriction like assumption 
(4.31) in Theorem 4.7 cannot be avoided. 

Example 4.11. Consider the case (4.34) with the function k(x) = § 
and xq = 1. According to (4.37), the trajectory of the Markov chain is de- 
terministic, with 



so that there is no explosion. This example illustrates that assumption (4.30) 
in Theorem 4.7 cannot be replaced by some arbitrarily slow growth at zero. 

4.3. Comments. In the pure coagulation case there is no explosion in the 
direct simulation model. According to Theorem 4.4, there is explosion in the 
mass flow model for a rather wide class of gelling coagulation kernels. 

In the pure fragmentation case there is explosion in both models. Theo- 
rems 4.3 and 4.7 cover wide classes of unbounded (at zero) fragmentation 
kernels. The sufficient conditions for the mass flow model are stronger than 
those for the direct simulation model. Example 4.10 shows that there are 
fragmentation kernels for which the direct simulation model explodes almost 
surely, while the mass flow model does not. Indeed, choose the kernel (4.34) 
with F(x) = x~ a ,a > 0, and k given in (4.39). Then the direct simulation 
process explodes almost surely for any initial value xq, according to Theo- 
rem 4.3. The mass flow process explodes with probability 1 — for any 
xq > i, according to (4.40). 

In the direct simulation fragmentation model explosion implies that the 
number of particles in the system reaches infinity in finite time (each jump 
increases this number). Due to mass conservation, this means that (at the 
explosion time) infinitely many particles are below any given size e, creating 
"dust." In the mass flow fragmentation model explosion occurs due to a fast 
approach to zero of a single particle. This explains what is going on in the 
example described above. In the direct simulation case each jump creates two 
fragments and some of them reach the interval (0, ^] leading to explosion. 



(4.40) 




For the choice F{x) 



x k = 2~ k , k = 0,l,.... 
logx + 1, one obtains [cf. (4.38)] 
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Fragments staying bigger than | do not avoid the fact of explosion. In the 
mass flow model there is only one particle, which stays above | with positive 
probability. 

Recently the interest in fragmentation processes has considerably increased 
(cf., e.g., [4, 15]). In this context the occurrence of small particles (dust) has 
been studied using different models. 

5. Concluding remarks. In this paper explosion criteria for jump pro- 
cesses with an arbitrary locally compact separable metric state space were 
established. These results are of independent interest. As an illustration, 
the general criteria were applied to stochastic coagulation-fragmentation 
models. The corresponding results (Theorems 4.3, 4.4 and 4.7) cover a wide 
range of coagulation and fragmentation kernels, for which phase transitions 
to infinitely large (gel) or infinitely small (dust) particles are known in the 
context of macroscopic equations. The proofs of these theorems are rather 
short, which illustrates the efficiency of the general explosion criteria. In 
particular, the results of Theorems 4.3 and 4.4 would hardly be available 
on the basis of the previously known criteria for one-dimensional processes. 
It might be of interest to extend the explosion results, which were obtained 
either for pure coagulation or for pure fragmentation, to more general situ- 
ations combining both processes and including source and efflux terms. For 
this purpose, perhaps more sophisticated choices of the test function r/ in 
the criteria will be required. 

A challenging problem for future research is the study of the limiting 
behavior (n — ► oo) of the stochastic coagulation-fragmentation models with 
explosion. In particular, the second (more difficult) part of the conjecture 
mentioned in Section 1 is still open. Convergence of a truncated mass flow 
model (particles exceeding a certain level are removed from the system) 
to the solution of the Smoluchowski equation has been studied in [9]. In 
general, the continuation of jump processes after the explosion point is a 
delicate problem. In the context of the coagulation-fragmentation models 
based on particles, there is a natural way of continuation — particles reaching 
or oo are simply removed from the system, while the others continue their 
evolution according to the previous rules. 

The conservation property of the mass flow equation [cf. (3.27) without 
source and efflux terms and for a constant test function] allows one to inter- 
pret its solution as a probability measure and to construct a related nonlinear 
Markov process. Such a process is determined by some stochastic equation 
with coefficients depending on the law of the solution. This approach has 
been carried out in [7] (pure coagulation), [17] (including discrete fragmenta- 
tion) and [13] (including continuous fragmentation). It would be of interest 
to clarify how the explosion phenomena observed in finite particle systems 
are represented in these limiting processes. 



EXPLOSION PHENOMENA 
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It is remarkable that some explosion phenomenon is probably recovered 
also in the direct simulation coagulation model in the limit n — > oo. A re- 
sult concerning explosion of an appropriately scaled tagged particle in the 
discrete case K(x,y) = (xy) a for 1/2 < a < 1 was announced in [21], but 
unfortunately has not been published so far. 
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